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Abstract 

A Comparison-Of. Two Typos of Velocity Models for the 
Lunar Crust: Smooth, Continuous and Stepwise Layered 

Anthony F. Gangi 
Department of Geophysics 
Texas ASM University 

The data from the Apollo- 14 and Apolio-16 Active Seismic Experiments 
have been reanalyzed and show that a power- law velocity variation With 
depth, v(z)»H0lV^ m/sec (0<z<10m) , is consistent with both the travel - 
times and amplitudes of. the first arrivals for source- to- geophone sepa- 
rations up to 32m, The-data were improved by removing Spurious glitches, 
filtering and stacking. While this improved the signal *tO-noise ratios, 
it Was not possible to measure the arrival times or amplitudes of the 
first arrivals beyond 32m. the data quality precludes a definitive dis^ 
tinctidn between the power- law velocity variation and the layered- Velocity 
model proposed previously. However, the physical evidence that the shallow 
lunar regolith is made up of fine particles adds weight to the 1/6-power 
velocity model because this is the variation predicted theoretically for 
self- compacting spheres. 

The 1/6-power law predicts the traveltime, t(x), varies with separation, 
x, as t(x) ■ t 0 (x/x 0 ) 5 / 6 and, using a first-order theory, the amplitude, 
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A(x). varies as A(x) - ^/ U ' m)/U . *U the layer-velocity model 
piediits t£c) * t 0 Cx/x 0 ) and A(x) * A 0 (*/x 0 ) 2 , respectively. The Measured 
exponents for. the arrival times were between 0*61 and 0.84 wliilc these for 
the amplitudes were between -1.5 and -2.2. The large Variability in the 
amplitude exponent is due, in part, to the coarseness with which the 
amplitudes are measured (only five bits are used per amplitude measurement) 
and the variability in geophone sensitivity and thumper-shot strengths. 

A least -squares analysis was devised which uses redundancy in the 
amplitude data to extract the geophone sensitivities, shot strengths and 
amplitude exponent. The method was used on the Apollo- 16 ASE data and it 
indicates there may be as much as 30 to 40% variation in geophone sensitivities 
(due to siting and coupling effects) and 15 to 201 variability in the thumper- 
shot strengths. However, because of the low signal-to-noise ratios in the 
data, there is not sufficient accuracy or redundancy in the data to allow 
high confidence in these results. 
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intwdneticri. the first liirar seismograms recorded by the Apollo- U- 
seismoueters (Latham, et at. , 1970a,b) surprised many seismologists. 

The ir-ttiusuol ly. long durations (see Figure U gave rise to msnerous 
theoretical speculatitms. Proposed mechanisms tanged from secondary- 
ejecta effects (Utham. et al.. 100a; Chang. et a!., 1970; HddmmedshauoV. 
1970) to scattering of the waves by shallow internal fractures and 
inhomogeneitics (Utham, et al., 1970a.b) or by topographic itregularrtres 
(Gold and Soter, 1970). It soon became clear that the secondary-ejecta 
mechanisms were not viable ones because the same long duration occurred 
for seismograms from amends with foci in the lmar interior. 

Early data indicated that the compressional-wave velocity was very 
low near the lunar surface (-0.1 km/sec; Latham, et al., 1970c, Sut 
and-Duennebier, 1970) and increased to approximately 6 Wsec at a 
depth. of 20 kilometers (Latham, et al., 1970d). Latham, et al. ( a, 
showe d that the variation of the aaplittaie envelope with time and distance 
was consistent with a diffusive-scattering mechanism provided the Q of 
the mediiift was greater than 3000. 

Gold and Soter (1970) interpreted the Apo, lo-l2 data to imply that 
the shallow lunar crust consisted of a deep layer of powder. They 
assumed a linear velocity Vatiatim With depth and, through computer 
simulation using ray acoustics, they were able to approximate the 
actual signal very well. They Showed that the long duration could be 
explained by scattering of the nearly vertically- incident waves by 
topographic irregularities (Figure 2). They also showed that the 
seismic amplitudes are greatly enhanced in such a mediua, so that it 
required less power to transmit seismic waves than previously believed. 


Kovach, et al- (l9n) proposed a layered model with a_stepwise- 
increas ing_yolocity variation, based on the. data of the Active Seismid— 
Eiqjeriment (ASH) at the Apolio-14 landing site. They obtained- a 
p & wave velocity (V p ) of 104 m/sec for. a top. layer Of 8.5 meters 
thickness and a V (|) « 229 m/sec for an underlying layer (the Fra Mauro 
formation) of 38 to 76 meters thickness. A similar model was used to 
interpret the Apollo-i6 ASE data and gave a V p * 114 m/sec for a 12.2- 
meter- thick top layer and a V p « 250 m/sec for an underlying layer 
TQ-meters thick (Kovach, et al., 1972). 

Gangi (1972) proposed a self-compacting-powder model which gives 
a velocity varying as the sixth root of the depth; in this model the 
velocity at the lunar surface goes to zero. This, in turn, gives a 
long duration to the signal by scattering from topographic irregularities, 
very low correlation between horizontal and vertical displacements , 
a changing signal envelope that varies with source- to- receiver separation 
and a varying spectrum over the signal duration. These effects have 
been noted by Latham et al. (1970c, 1970d) and they also are explained 
by the diffusive scattering model (Latham, et ah,1970c) and the surface- 
irregularity scattering model (Gold and Soter, 1970). 

Kovach and Watkins (1973) extended and refined the layered model 
by incorporating the traveltime of the Apollo- 14 Lunar-Module ascent. 
However, they pointed Out that; "the exact details of the Velocity 
variation in the upper 5-10 km of the Moon cannot yet be resolved 
(i.e., whether it is smooth as depicted or a stepwise increase) but one 
simple observation can be made. Self- compression of any rock powder 
such as the Apollo 11 or 12 soils or terrestrial sands cannot duplicate 
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the observed magnitude, of the luiar velocity change and the steep 
velocity-depth gradient (-2 km .‘W However, .it U-not expected 
that a self-COmpacting-powder layer of 5 km thick Would.exist oh the 
Moonu if such a layet exists, it would be, most likely, thinner than 

1 tort arid probably thinner than. 100 nu 

Dainty, et al. (1974) performed a detailed analysis of the diffusive- 

scattering mechanism and compared their theoretical results both with 
lunar data and seismic-modelling data. They showed they could match 
the envelopes of the luiar seismograms using this theory if, for a 
frequency of 0.4S Hz, the apparent thickness of the scattering layer 
is 25 km, the mean distance between scatterers at the base of the layer 
is -5 km and the Q of the medium is 5000. The corresponding values 
for a frequency of 1,0 Hz arei 14-km scattering- layer thickness, -2 km 
between scdtterets and a Q of 5000. The thicknesses of the scattering 
layer (and its variability with frequency) seem to be inordinately large 
and may indicate that the model used is not appropriate for the lunar 
crust. A similar analysis should hold for body-wave scattering by 
topographic irregularities-, in this case, the scattering- layer thickness 
would correspond to the surface area over which the nearly vertically- 
iiiclden t waves are efficiently scattered and the spacing between 
scatterers in the layer would correspond to the spacing between surface 

scatterers (of Wave-length size). 

Cooper, et ai. (1974) used the data of the Active Seismic Experi- 
ments of Apollo 14 and 16 along with the Lunar Seismic Profiling 
Experiment (LSPE) data 6* Apollo 17 and other man-made impacts to 
obtain a model of the velocity structure of the shallow lunar crust. 
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They assumed a layered model and assumed that the first arrivals- 
(beyond about Id m) were seismic refractions. They found- their travel - 
time data, were consistent with a f iveJayer model in which the. velocity 
is: 1) 100 m/sec in the top. layer Of 4-m thickness - t 2) 327 m/sec in 
the. next layer to a depth of 32 m (thickness Of 28 ui)* 3) 495 m/sec to 
a depth of 390 m, 4) 960 m/sec to a depth Of 1385 m and 5) 4700 m/sec 
for a depth down to at least 1800 m. However, this last velocity_iS 
determined from a single source (the IM impact) at distances of the 
order of 8.7 km from the geophone array (four geophones) . The shallower 
structure is obtained from the traveltime data resulting from the 
eight explosive-package detonations and the LM ascent; all these sources 
are within 3 km of the geophone array* Cooper, et al* (1974) show these 
data can be fitted well with a continuous, linearly- increasing Velocity 
with depth, z; namely, V ** 395 + 778 z (m/sec) for z in meters. They 
also state that "Various power law velocity models can be made to fit 
the observed data ..." when only the explosive -package and LM- ascent 
data are Used. 

it is clear there is still Some question regarding the velocity 
Variation with depth in the shallow lunar crust (z<l km). Since the 
shallow lunar crust severely modifies the received Signals, evert those 
from large distances, it is important to know this shallow velocity 
variation well. Therefore * it is worthwhile to reanalyze the data to 
determine which velocity variation with depth is the most probable. 

The data from the Apollo-14 and Apollo- 16 ASE’s have been reanalyzed 
the the results ate given below. 
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Apollo- 14 and -16 ASB Data. The data used in this analysis- are from 
the astronaut-activated thumper device, of thcJ\pollo-14 and,Agpllo*16 
ASE's. - In both experiments * three geophones were sited on the surface 
in a linear array with 45*72 m (150 ft) spacing_between geophones 
(Lauderdale and Eichebnan, 1974). The thumper device was fired at. 

4* 57 m-(15 ft) intervals between the ends of the arrays (see Figure 3) . 
Firings (shots) 5, 6,.. 8, 9, 10, 14, 15 and-16 of the Apollo- 14. ASE 
misfired and.no data are available for them* For the Apollo- 16 ASE, 
two shots were omitted between -geophones 1 and 2; namely, those at the 
4.57 m spacing from the two geophones. 

The signals from the geophanes are sampled every 1*887 msecs, 
corresponding to a Nyquist frequency of about 265 Hz. Because of data 
transmission limitations* a trade-off between sampling rate and the 
number of bits per sample had to be made. The result was that only five 
bits were available for each sample. In order to cover the maximum 
possible dynamic range with only 32 possible binary numbers* the seismic 
signals were log Compressed for large signal levels. The correspondence 
of the binary-data values (0-31) and the voltage from geophone 1, Apollo 
16 is shown in Table I. TTie other geophone voltages have similar 
correspondences with the binary data. With only 32 levels possible 
for the geophone output voltage* the resulting traces will have a coarse 
character. This makes it difficult to obtain accurate amplitude informa- 
tion if no processing or filtering is performed On the data. Fortunately, 
it is possible to process the data to obtain reasonably accurate 
amplitude values. 


In ordet to- achieve meaniful results from the analysis, it was 
necessary to improve the otiginaUSE data. Figure 4 shows three, 
representative traces of the taw data from the Apollo- 16 A$E.- These 
data-are from the tenth thumper shot and the soutCe-to-reCelver 


separations are 50.2® m (165 ft), 4.57 m (IS ft) and 41.14 at (155 ft) 
for geophones 1, 2 and 5 respectively. The thumper- firing time is 
1.2 seconds after the beginning of the traces. While a high signal - 
to-noise ratio (S/N) exists for the shortest separation, the S/N for 
the other two traces is so low that it is difficult, if not possible, 
to pick the first arrivals or to measure their amplitudes. In addition, 
geophone 1 shows severe "glitches", most of which are almost uniformly 
spaced in time and of uniform anpHtude but there are others of varying 


amplitude and times of occurrence. Similar large glitches are seen on 
the other two traces. A close look at the data showed that there are 


cellar glitches throughout the records; these are recognized by the 
fact that they are of short duration — generally, only otte or two 
samples - and had values which were inconsistent with proceeding and 


following ssurtpie values. 

The first data-improving operation performed was to go through 
the data by. hand and remove the extraneous values and replace them by 
values interpolated from neighboring values. A computer program was 
not used in this process because; 1) there are relatively few glitches 


(excluding the regular, perodic ones in geophone 1, there are fewer 
than It), 2) the coarseness of the amplitude values precludes automatic, 


computer interpolation and 3) a mmter of different criteria wore 
used simultaneiously to identify and correct the bad Sample values. 
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The result of the "deglitching” process is shown, in Figure 5 for. the 
same traces shown in J&gure 4. While this improved the records cent 
s idcrably , it is clear the S/N ratios for the gcophone- 1 and -3 traces 
are still-too low to allow positive identification of the first arrivals. 

TO improve iha S/N and smooth out the traces, the data were bandpass 
filtered with a four-pole , anti-aliased, Butterworth filter (-12 db/ 
octave slopes at both iow and high frequencies) which had 3-db fre- 
quencies at 10.5 Hz and 66.25 Hz. The result * for the same three traces, 
are shown in Figure 6. While this improved the S/N significantly and 
improved the character of the traces (compare Figures 5 and 6), the 
S/N for separations larger than 9.14 m (30 it) was still low because 
of the decrease in the direct arrival’s amplitude. 

Spectral analyses of the seismic traces were made to determine the 
frequency bond of the seismic energy and to see if there was significant 
aliasing of the data. The amplitude spectrum (for geophone 2, shot 10* 
Apollo- 16 ASE) of the first two seconds (1024 samples) is shown in 
Figure 7. While only half of the full amplitude spectrum (0 to 265 Hz) 
is shown there, it is clear that there is little, if any, aliasing 
because most of the signal energy is contained between 10 and 90 Hz 
with the major part between 10 and 40 Hz. This is the spectrum of the 
middle trace shown in Figures 4, 5 and 6. 

To further improve the data* the traces with the same source -to 
receiver separation for both ASE' s were summed (or "Stacked") together. 
The implicit assumptions being made here are: (1) the velocity 

variation with depth is the same at both the Apollo- 14 and -16 sites 
and (2) there is lateral homogeneity for the direct waves at both 


sites. The first assumption is reasonably consistent with the results 
found by Kovach, et ai. (1971) for the two- sites', the second assumption 
is consistent with the. equivalent assumption made by Kovach, et ai. 

(ld?i) in their interpretation Of the data at each site,- 

The traces that had the same source -to- receiver separation are 
listed in Table II for both the Apollo-14 and -16 ASE’s. The thumper- 
shot numbers, corresponding to the given shot -to- geophone separations, 
are listed in the right half of the table. Among the two experiments, 
there were between 4 and 7 traces with the same separation. If the 
background noise is random and the assumptions cited hold, the stacking 
should give improvements in S/N between 2 and /7. The resulting sum 
signals were amplified so that the peak excursions would be plotted 
almost full scale for each trace. A representative result is shown 
in Figure 8. The second trace in Figure 8 is at the same separation 
as the middle traces in Figures 4, 5 and 6 (i.e., the geophone- 2 
trace for the 10th thumper Shot Of the Apollo- 16 ASE). For this trace, 
the S/N improvement should be better than a factor of 2; however, this 
degree of improvement was not achieved. Nevertheless, improvements 
‘n S/N were achieved for this trace, and for the other traces at larger 

separations, by the stacking technique. 

The result of summing the deglitched traces is shown in Figure 9. 
These signals were filtered, before summing, with a 4-pole, anti-aliased, 
bandpass, Butterworth filter with 3-db frequencies at 20 and 50 Hz. 
Arrival times can be determined with some certainty for sopaiations 
up to 32.00 m (105 ft); it is difficult, if not impossible, to pick 
arrivals beyond that distance. 
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One o£ the single-geophorte profiles (geophone 1, Apollo 16) is 
shown in Figure 10. Arrival times can be easily picked for separations 
up to 18.29 m (60 ft) and, with difficulty, for 22.86 m (75 ft) and 
27.43 m (90 ft). At 32.00 m separation (105 ft), the first arrival 
is burl* in the noise. We were net able to determine a first arrival 
with any degree of certainty for separations greater than 32.00 m 
(105 ft). This is consistent with the finding of Kovach and Watkins 

(1973) for the thumper shots. 


Results. The traveltimes and amplitudes of the direct (first arriving) 
seismic signals of the Apollo-14 and -16 ASE's were analysed. The 
traveltimes and amplitudes for separations up to 32.00 m (10S ft) were 
obtained both from the Surted (stacked) traces and from individual 
traces. In one case, all the "noise-free" traces from both ASE's were 
stacked in on attempt to improve the S/N ratio. In two other cases, 
only the "noise-free" traces from each ASE were stacked to give Apollo 
14-only and .Apollo- 16-only stacked profiles. If there are significant 
differences in the velocity structure at the two sites, these individual- 
site stacks would show the difference. Little difference was found, 
over the 32 m, in the traveltimes for these two stacks. The traveltimes 


for individual geophone profiles Were also measured to test the ass«p- 
tioh of lateral homogeneity at each site. Ihe quality of the data 
precluded any positive conclusion regarding this assumption; however, 
the improvement in S/N ratio achieved by the various stacking indicate 
this is a reasonable assumption. 
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Tray alt imes . The traveltimes- for five of the cases investigated ate 
listed in Table II L In these cases whete the S/N ratio was high (up 
to and including 18.29 m separation*, the traveltimes could be determined 
to within 1/2 sample time (+ 1 msec) . However,, systematic errors — 
such as these due to variations in the separations* elevation differences, 
shot- times* etc. — could be as high as one or two sample times. 

Log/log plots of traveitime versus separation were used to test 
the hypothesis of a power- law velocity variation. It can be shown, 
using Kaufman's (1953) work* that a velocity variation with depth, z, 
given by 

v(z) = v 0 (z/z Q ) n (1) 

results in a direct -wave traveltime, t(x)* with separation, X* given 
by (S6e also, Gangi, 1972) 

t(x) • t 0 (x/x o ) 1 ' n (2) 

where t Q is the traveltime corresponding to the separation x Q and v Q 
is the velocity at depth z Q . This incorporates both the traveltime/ 
separation relationships for a constant-velocity medium (n*0) and that 
for a Self-compacting-pOwder medium (n=l/6) . Therefore* for the two 
power- law velocity models (n=0 and n=l/6) , the traveltime curve in a 
log/log plot would be a straight line whose slope* m, would be determined 
by the power- law exponent (m=i-rt). 

The slopes of least ‘Squares -fitted straight lines are given in 
Table III along with the velocity v Q which corresponds to the velocity 
extrapolated to z Q “l km. As indicated earlier, it is not expected that 
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the powdered layer would extend to 1 ton; therefore, v 0 is not an 
estimate of the velocity at that depth but is merely a constant used 
to characterize the .velocity. IV depth z 0 *l toil is chosen only for 
convenience; the reference depth could have been chosen to be 1 m, 

1/6 

iii which Case, the v 0 's in Table III would be multiplied by (.001) 

* 0.3162. While the measured slopes are variable, they are all 
consistently lower than mFl-n=i, the value that would be obtained 
for the constant velocity model. Hie measured values tend to cluster 
near the Value predicted by the self- compact ing - powder model; namely, 
111*1-1/6=. 833. 

The variation in the reference velocity, v Q , is much greater than 
that of the Slopes; its values vary between 340 and 630 m/sec. The 
slope of 0.63 and reference velocity of 1200 m/sec for the Apollo-16, 
geophone- 1 profile (column 5, Table III) are not very accurate because 
there are only three good data points (the traveltimes at 9.14, 13.71 
and 18.29 m) for determining these values. It gave the least consistent 
values for n and v 0 . in computing the least-squares lines, the question- 
able data Were given a weight equal to one-quarter that of the high-S/N 
data. 

Traveltimes were calculated from the Apollo 5 - 14 and Apollo- 16 Velocity 
models given by Kovach and Watkins (1973). These are shown in columns 
B and C of Table Hi. in column A, the traveltimes for a powder* layer 
model with v o *350 m/SeC and n*l/6 are tabulated. This latter model was 
an average model found from all the cases treated when the velocity 
exponent, n, was constrained to be 1/6. Overall, there are not large 
differences between the measured traveltimes and the calculated travel- 
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tiroes using any o t the models. However, the biggest differences between 
the Kovach end Watkins-models .and the measured. Values occut-at the small- 
separations, .precisely where the $/tL ratios are highest. and where the 
traveltimes can be picked with the greatest certainty. Their models 
can be made to fit the close-in data simply by introducing a thin> lower 
velocity layer at the surface. But it should be. recalled that they 
already have low velocities for the top layers (104 and 114 m/sec fot 
Apollo- 14 and -16, respectively) which are relatively thin (8.5 and 

12.2 m, respectively). 

The traveltime data for the combined Apollo- 14 and -16 stacked 
traces (column 1, Table III) are shown in Figure 11 along with the 
least -squares -fitted line. These data are from the deglitehed traces 
which have been bandpass filtered with a fourth-order, Ruttcrworth 
filter having 3-db points at. 3 and 66 Hz~_ It can be seen fronuFigure 
11 that the straight line is an excellent fit to. the data and that it 
would be difficult to change the slope from its given value (0.76) 
to 1.0, the latter value corresponding to the constant -velocity model. 
Equally good fits of data points to straight lines were found for the 
Apollo^ 14 -only and Apollo -16 -Only stacked data. 
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Am plitudes , llie traveltimes of the-first arrivals over the 0-32 m range 
do not demonstrate a clear distinction, between the powered- layer and the 
layered-velocity models. The data. accuracy is such that cither-model can 
be accepted* . To try to distinguish between the two models, the amplitudes 
of the first arrivals were measured and compared with the expected distance 
variation predicted by the two models. 

Since the thumper shots give primarily vertical forces and the geophones 
are vertically oriented, the amplitude of the direct p-wave arrival in the 
layered model should vary as the inverse square of the separation, 

AW - A 0 (x/x 0 )- 2 4 { j, 

for small separations (see, for example* White, 1965, p. 215). On the other 

hand, for a power- law velocity model, the amplitude variation with separation 
is given by (see Appendix A) 

( 4 ) 

where S(p)dp is the energy radiated in a bundle of rays having ray parameters 
lying between p-dp/2 and p+dp/2, the ray parameter is given by p»sin0 (z)/v(z), 
6(z) is the angle between the ray and the vertical (z) direction, v(z) is the 
velocity variation With depth and t is the traveltime for the ray (with ray 
parameter, p) which returns to the surface at separation x. For the self- 

compact ing-powder model* the amplitude variation is estimated to be (see 
Appendix A) 


AW - A o (x/x b )-9 3 -”0/12 i m>1> 


( 5 ) 


where in is a measure of. the source radiation-pattern in the^ power- law* velocity 
medium, to insure integrability of 

E - $ S(p)dp, 

where E is the energy radiated by the thumper source, we find m>l (see Appendix 
A). This indicates the amplitude decrease of- the direct wave with separation 
is less in the powder- layer model then that in the constant- velocity model. 

This is consistent with, the conclusion of Gold and Soter (1970) based on 
their analysis for a linearly increasing Velocity with depth. 

The determination of the amplitude variation with separation for the 
Apollo- 14 and -16 ASE data is more difficult than determining the traveltime 
data because of: 1) the coarseness of the amplitude sampling * 2) the 

variability of the thumper-shot strengths, 3) the variability of the geOphone 
se nsi tivities (primarily due to siting and coupling, of the geophones) and 
4) the low S/N ratio for the larger separations. The- coarseness of the 
amplitude data is significantly reduced by the interpolating effect Of 
bandpass filtering. The variability due to the shot strengths, the geophOne 
sensitivities and the low S/N ratio are reduced by the averaging inherent 
in stacking or surrinihg traces (provided the signals are sufficiently coherent 
for a given source/receiver separation). 

On the basis o* the measured arrival times (at least for separations 
less than 22.86 m - see Table III), sufficient coherency Of the signals exists 
s. > that averaging of the amplitudes should be possible by summing of traces. 
The measured amplitudes are given in Table IV. Both the amplitudes for 
individual truces and for stacked truces are given. Measurements were made 
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on data that had-been bandpass filtered by anti-aiiased, .fourthiorder - 
Butterworth filters with 3 db frequencies of. 3 to 66 Hz and 20 to 40 Hz- . 

It can be seen that there is a great deal of scatter, in the data_ Some of 
this is diie to the thumper- shot variability and some -due to geoghone siting, 
but the major part is due to low S/N ratio and the coarseness Of the 
amplitude data. Straight lines were fitted, by least squares, through the 
data points (an a log- log graph) and the slopes of these lines are included 
in Table IV. A representative plot of the amplitude data along with its 
least-squares - fitted line is shown in Figure 12. This represents one of the 
most complete sets of amplitude data available for a single geophone; namely, 
geophone 3 for the Apollo-16 ASE. The original traces were bandpass filtered 
(3 db frequencies at 3 and 66 Hz) prior to measuring the amplitudes. 

Because of the low S/N ratio at the larger separations, it is not 
certain that a straight line (oh a log- log plot) is the appropriate fitting 
function. While all the data are fairly well fitted, by the line in Figure 12 
(with a slope equal to -2.01), it is clear that the two largest amplitude 
values (at 4.57 and 9.14 m), which have the best S/N ratios, suggest a lower 
slope. 

The slopes found for all the cases with fairly good data lie between 
-1.5 and -2.1. However, the possible errors on these slopes are of the order 
of ±0.5. The fact that the slopes are more negative than -1 and close to -2, 
the slope predicted by a simple flat- layer model, does not mean the amplitude 

tgjf 

data verifies that model. From equation (5) , the slope predicted/the powder- 
layer model would be more positive than -1. However, this equation and the 
theory used to predict a slope of -2 for the flat- layer model are based on 
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siaplyfying assumption; namely, that all the-sources are of equal strength, 
the geophones are equally coupled to the regolith. therejs no attenuation by 
absorption in either model, there iS-no energy loss byconversion of p-wave 
energy into-s-«av«nergy(for the powder-layer model) and there -are no 
scattered in the lunar regolith. The latter three effects would increase 
the -amplitude loss with distance so that the predicted slopes (-2 for the 
flat layer and - (13-m)/l2 for the powder layer) should be considered upper 
bounds on the measured ones. The variability of the thumper-shot strengths 
and of the geophones would increase the scatter in the data. 

While the aaplitude data do not preclude either model Conclusively (as 
they would have if the measured amplitudes decreased more slowly than inversely 
with separation), they do favor the. powder-layer model. All the loss 
mechanisms lead to a greater decrease in aaplitude than_p.redicted by the 
simile analyses of the two models. However, the aaplitude data do not show 
a more rapid decrease than that predicted for the homogeneous-layer models 
proposed by Kovach and Watkins (1973) while the data clearly do show, a more 
rapid decrease with separation than that predicted by the simple (first-order) 
theory for the powder-layer model. This discrepancy in the amplitude 
variation with distance can not be explained by interference of other Waves 
with the direct wave. For the short separations where aaplitude data is 
available (generally less than 27.43 m), interference from reflected or 
refracted waves would not affect the amplitudes by interference for the flat- 
layer models; nor would a Velocity discontinuity at a depth greater than 

about 10 m effect the aaplitude results (by the Same types of interference) 
in the powder- layer model. 


Gconhono-Counl ing and Shot-Strength Variability , to eliminate the effects 
OC-variability in the. geophone coupling. and the thumper “Shot strengths, an 
analysis of the amplitude data was made which determines both the geophone 
sensitivity (in-place) and the. shot strengths as well as the exponent of 
the amplitude variation when there is sufficient redundancy in the data. 

If it is assumed that either the flat- layer model or the powder- layer 
model is valid, the measured amplitude at a particular geophone due to a 
particular source will be given by 

C6) 

where is the sensitivity of the i-th geophone (including coupling and 
siting effects) iocated at x i , is the strength of the j-th shot located 
at and m is the. exponent of the amplitude variation. 

Equation (6) Can be normalized to the sensitivity Of a particular 
geophone, say Gj (1=1,2 or 3), and to the strength of a particular shot, 
say Sj. This normalization is necessary because, quite clear iy, each 
geophone sensitivity can be multiplied by some constant factor and each shot 
strength divided by the same factor without changing the resulting amplitude. 

Letting GjSj“A 0 i equation (6) becomes 

Ay - A 0 (G 1 /Gj) (Sj/SjJlxj-Xjf. (7) 

Equation (7) can be linearized in terms of the relative geophoiie sensitivities 
the relative shot strengths, the exponent m and the arbitrary constant by 
taking its logarithm: 
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log Ajj ■ log A d + log(G |/Gj-) + logCSj/Sj) + rtlog(x^*x^| 
oiTj for convenience in writing, 

a ij “ a 6 + h * s j * ^ij 


( 8 ) 


vAiere a.^iog A^, X^loglx^x. | , g^logCG./Gj) and sflog{S./$j). 

The optimum values, in a least-squares sense, of a^ g^ and m can 
be determined by minimizing the summed, weighted and squared error 

2 - 1 J 

E (V**’ 5 ) =J i ^ w ij f a 0 + % + «i + s j - a i;j ) 2 , (9) 

as a function of these parameters. The result (see Appendix B) is the 
matrix equation 


a * A • p 


(io) 


where a is a vector whose components depend only upon the measured 
amplitudes the weights (w^j) and the measured separations (X^ ) ; 

A is a square matrix whose components depend only upon the weights and the 
measured separations while p is a vector whose components are the unknown 
parameters: g^..., s^..., a^ and m. (The detailed form of this equation 
is given in Appendix B). The solution to this matrix equation is 


_ *-i ^ 

p * A *a. 


(ID 


If there is sufficient redundancy in the data, the matrix will be well 
conditioned and non-singular and will have a stable inverse. 

The weights are established from the quality of the data. The Weights 
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for. the Apollo- 16 data are shown, in table V. From the table it is seen- 
that only 14 of the 19 thumper shots gave useful amplitude data (shots 1, 

U, 13*. 14 and. 19 wore not useable) and, of those 14* only throe (shots 6,_ 

7- and 17) give measurable first-arrival amplitudes on more than one geophone. 
(Shots 12 and 17 gave amplitudes of .89-and..06?, respectively, for the 
3 s - 66 Hz bandpassed traces , on geophone 2; all other amplitudes are given in 
Table IV). Therefore, only, six-amplitude measurements (two for each shot) 
are available to determine the six parameters a Q , m, gp gp s 6 , s.^ (when 
geophone 2 and shot 1 are used as the reference geophone and shot, respectively). 
With the relative geophone sensitivities, the constant a Q and the exponent 
m set by these data, the remaining relative shot strengths will be determined 
by the assumed amplitude variation (equation 6 or 7) and the measured amplitude. 

Having only six correlative amplitude measurements to determine six 
unknowns (by means. of the linear equations 6 or 10) means there-is little 
redundancy in the amplitude data* Nevertheless, the solution of these six 
equations in the six unknowns do constitute a ieast-squares solution. This 
is because weighting factors are used in the equations; the weights can be 
interpreted to mean that more than six measurements of equal weight were made, 
some of which were identical measurements (i.e. , Same shot location), and 
the results combined together to give a single result of greater weight. 

Using the six available correlative amplitude values, the matrix 
equation becomes: 
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<w ij 4 ij > ' 



V 

u. 

V 

<w 3J > 

<w i6 > 

<w i47 > " 

*0 



m 

•iWijX i} > <w tj x i;j > 

<w 3j X 3J > 

<w i6 X u > 


it i 

<W lj a ij > 


m 

• “V 

0 

<w 16 > 

^1,17* 

8 1 



- 


< V 

<w 36 > 

1 

<w 3,17 !> 

> 

8 3 

<W i6 a i6 > 


1 

“ m 

- 

<w i6 > 

0 

s 6 

<w il7 a il7 > 



- 

- 

- 

<*1,17* 

S 17 


where i 1, 2 and 3 and j B 6, 7, 17. In terms of assumed values 
of w.. and the measured values of and a.^ this equation becomes 


-3.734 


3.50 

10.35 

i.oo 

.75 

1.00 

1.25 

-H.710 


- 

30.84 

2.62 

2.39 

3.13 

3.48 

-0.211 


- 

- 

1.00 

0 

0 

1.00 

-1.164 


- 

- 

- 

0.75 

0.50 

0 

-1.592 


- 

- 

- 

- 

1.00 

0 

-0.914 


- 

- 

- 

. 


1.25 


*1 

s 6 
. s 17 


Solving this matrix equation, the relative geophone sensitivities 
and relative shot strengths are found to be: 

* .724; “1.40 

V S 7 * *803; S 17 /S ? - .848 


<iiifc¥irwntfrrif jjjfc w 1 JrtnBBwaa rni 1 1 ~ 1 1 
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and the exponent is 

m * -3.57. — 

Unfortunately, these values appear to be unreasonable; this is not 
surprising considering the lack of redundancy and quality in the amplitude 
data. The 30 to 40% differences in the relative geophone sensitivities are 
not too unreasonable* but ate higher than expected. Also, the 15 to 20% 
variations in the shot strengths are possible, but again seem large. The 
value of the exponent (hf-3.57) is different by almost a factor of two 
compared to the values obtained using single- geophone profiles and stacked 
profiles (compare Table IV). The 30 to 40% differences in geophone 
sensitivity have no effect on the amplitude Variation with distance as 
determined by a single-geophone profile. The 20% differences in shot 
strengths (of shots 6 and 17 relative to shot 7) would not Cause appreciable 
differences in the slopes (or exponents) obtained from single-geophone 
profiles (provided, of course , that these differences are representative 
of the differences in the other shots) . It is concluded that the least- 
squares analysis given above does not give reliable values for the parameters 
(m, G 3 /G 2 , S 6 /S 7 , S 17 /S ? ). However, the method is a valid one and 

the reason for the unreliability in the parameter values is the lack of 

redundancy and quality in the data. 

While the method is not useful for this data Set, it is presented in 
detail because there may be other instances where it would give valid results. 
It provides a rationale for the design of seismic experiments which test 
amplitude variation with separation when variability in source strengths and 
geophOne sensitivities is anticipated (as is generally the case). 
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The same amplitude analysis could not be performed on the Apollo- 14 
ASE data because there were no correlative amplitude values for geophones 
2 and 3 (due to misfires and poor signal-to-noise ratios) . 
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Summary 

The data from the thumper shots of the Apollo- 14 ahd Apollo- 16 ASE's 
have been reanalyzed to test whether the velocity variation in the shallow 
lunar crust (depths < 10 meters) can be represented by a self-compacting 
powder layer as proposed by Gold and Sotor (1972) and Gangi (1972) or by 
constant velocity layers as proposed by Kovach et.al. (1971, 1972, 1973). 

Both the traVeltimes and the amplitudes for the first arrivals were 
remeasured and concaved with the values predicted by the self-compacting- 
powder- layer model proposed by Gang! (v(z)®v 0 (z/z 0 )^^) and the layered- 
velocity model proposed by Kovach et.al. To improve the quality of the 
data, they were "deglitched" to remove spurious values and bandpass fil- 
tered. Four-pole, anti-aliased Butterworth filters with bandpasses between 
3 and 66Hz and 20 and SOHz (3 db frequencies) were used to improve the 
signal- to-noise ratio (S/N). In addition, traces from different thumper 
shots and with the same source-to-geophone Spacing were summed together to 
improve the S/N. While these techniques improved the S/N, it still was not 
possible to measure traveltimes or amplitudes of the first arrivals for sep- 
arations greater than 32m. 

While there is variability in the results obtained (see Table ill) , the 
traveltimes for the direct arrival over a separation of 32’m can be fit by the 
1/6-power velocity model. Hie measured values of the exponent for an assumed 
power- law velocity varied between approximately 1/3 to 1/7; that is, 

.67 <l-n <.86 (see Table 111) where n is the exponent for the depth variation 
of the velocity. Hi# best (or average) model for both the Apollo-14 attd Apollo-16 
sites is v( 2 )« 350 (z/z o ) 1 ^ 6 m/sec for z Q *l km or v(z)«110z 1//6 , 0<Z<10m. This is 



- -^ B ra n niw 




{.Jy. (! 
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fcirly close to the velocity variation, v (z ) * 190z 1/6 *_ predicted by Gattgi 
Cid72) on the basis of Gassmannls analysis (1953) and the measured mechanical 
properties of the lunar soil* _ 

The measured traveltimes of the first arrivals over the 32m. separation 
are in reasonable agreement with the values predicted by the layered model 
(see Table III). However , the biggest percentage deviations occcur at the 
two Shortest distances (4*57 and 9.14m) where the S/N is high . and the 
traveltimes can be measured most accurately* At these separations, the 
measured arrival times* which are accurate within at least one sample interval 
(or 1.89 msecs), differs from those predicted by the layered model by 10 to 
15 msec. The corresponding differences for the power-law model is generally 
less than 2 msecs. While this indicates that the self-compacting-powder- 
layer model is probably the correct one, the quality of the data precludes a 
definitive distinction between the two models. 

No comparison was made of the measured traveltimes with those predicted 
by the linear Velocity variation used by Gold and Soter (1970), namely, 
v( z )=v 0 +az, because it was an assumed velocity variation which is hot based 
on any physical mechanism. The traveltime relationship for this velocity 
variation* t=(2/a)sinh (ax/2v Q ), should also fit the data to the same accuracy 
as that of the layered- Velocity model. It, too, would have the largest 
percentage deviations at the shortest distances. 

Ah analysis of the amplitudes of the first arrivals was performed to 
test the models. Hie predicted amplitude variation with Separation, x — as- 
suming no amplitude loss due to attenuation, scattering or conversion of p-wave 
energy into s-wavc energy - for the layer model is x* 2 while that Tor the 
1/6-poWer Velocity model is m)/12^ The measured exponent varied 
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from -1.55 to *2.34 (see Table IV) with the average value near *2.0. While 
this result, at first glance, seems to -favor the Constant velocity model, 
the fact that there will be amplitude, less due to scattering, attenuation 
and wave* type conversion makes this result more consistent- with the power- 
law model. However, the large errors in the amplitude data - which are more 
severe than the errors in the arrival times - preclude a definitive 
conclusion regarding which is the appropriate velocity mode/. 

An attempt was made to eliminate the errors in amplitude, due to vari- 
ations in the geophonc sensitivities and shot strengths, by using a least- 
squares method. The method requires that the signals, from individual shots* 
be detected on two or more geophones. Unfortunately, only three thumper shots 
were detected On pairs of geophones, and no thumper shots gave detectable 
first arrivals on all three geophones. Consequently, there was too little 
redundancy in. the data to give reliable values for the relative, geophone 
sensitivities, relative shot strengths or the exponent for the amplitude 
variation with separation. Only for the Apollo- 16 ASF. was there sufficient 
data to perform this analysis at all, and it indicated that there could be 
30 to 40% variability in the geophohe Sensitivities and 15 to 20% variability 
in the thumper-shot strengths. Ah amplitude variation with Separation equal 
to x" 3 * 6 was obtained from this analysis. It is not possible to give much 
credence to these values because the amplitudes used in this analysis were 


small and had large variability. 

In conclusion, it has been demonstrated that the power- law- velocity 
model predicts: 1) the measured arrival times of the first arrivals as well 
as, if not better than, the layered-velocity model does and 2) the amplitude 


variation, with separation as-weil as that model does, the quality of the 
data does not allow* definitive choice, to be-made between the_two.models.- 
However, the power- law model predicts a very low velocity at the lunar six r f &ce 
which, in turn,, implies that seismic rays will be nearly normally incident to 
the surface. This would explain why there is little. correlation between the 
Vertical and horizontal components of the motions detected by the Passive 
Seismic Experiment seismometers. It also implies that the long duration of 
the seismic signals detected on the Moon is due to scattering by even shallow 
undulations of the surface (Gold and Soter, 1970 and Gaigi, 1972). The power- 
law velocity model also predicts that the lunar regolith is composed of fine 
particles (Soil) down to a depth of S to 6 meters. The power- law model 
indicates that the velocity below 6 meters is not 'Sampled" by the first 
arrivals detected over separations less than or equal to 32m. 
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Appendix A 

Variation in amplitude with distance in a vertically inhomogeneous medium— 

An approximate analysis of the variation in amplitude of a compress ional 
wave in a vertically inhomogeneous elastic medium can be made using ray 
theory. The analysis closely follows the developments given in Bullen (1963) 
arid Officer (1958). 

The analysis is approximate in that it does not take into account 
either the variation in Waveform of the propagating wave (i.e. , dispersion) 
or the conversion of p-wave energy into s-wave energy (these assumptions 
are also made in the above references). 

We assume that, for a source on the surface, the energy, dE, contained 
in a "bundle of rays" with ray parameters between p-dp/2 and p+dp/2 is equal 
to the intensity (or energy per unit area), I, times the area subtended by 
the ray parameters (see Figure A. la) 

dE(p) = I(x\p)dA = S(p)dp (A.l) 

where 

I (x* ,p) * the wave intensity for a ray with ray parameter, p, at a 
horizontal distance, x' , away from the source point, 
dA * area contained between the circular "cones" given by 
p-dp/2-COrtstant and p+dp/2«constartt (dA=2irx , dw) , 

$(p) ■ energy per unit change in the ray parameter, 
p ■ the ray parameter ■ sin0(z)/v(z) , 

0(z) ■ the angle between the ray and the vertical z axis and 
is measured counterclockwise from the z axis, and 


v(z) * the velocity variation with depth. 

The intensity for any point along the fay can be expressed as 


,p) 


- SCP) dr 


(A. t) 


At a fixed depth, z, and for a particular ray bundle centered about ray 
parameter p, (see Figure A, lb) 


dw dx‘ 

dw(x' ,p) - cos 0dx' ; or 3- ■ cos0 


(A* 3) 


(A. 5) 


When the ray reaches the surface (z=0), 0(O) a ir (cos0 (0)--l) and 

x' (piZ=0)=x(p). *^ A ’^ 

The intensity at the surface receiver then becomes 

Kx.p) * -ffe 5 $r. <A ‘ 5) 

(Note, the intensity is positive since both S(p) and x are positive but 

dp/dx is negative at the surface). 

The general relationship for vertically inhomogeneous media 

p - dt/dx (A - 6) 

holds, where t(x) is the traveltime; therefore we have 

. SfD) d 4 t (A. 7) 


_ f , S(p) d 2 t (A. 7 

= m dr. 

It can be verified directly that p-dt/dx for a velocity variation 


of the form 


v(x) = v 0 (z/z o ) n (0<rt<l) 
by using: (1) the fact that 

dt/dx * (dt/dp)/(dx/dp) ^ A ‘ 9 ^ 

and (2) the parametric equations for the traveltime, t(p), and the source/ 

receiver Separation, x(p), (see, for example, Kaufman, 1953) 


(A. 9) 
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t(p) & c n p" (1 ' a)/n 
x(p)_* (l*n)C^p' 1/n 
Where is a constant equal to 


*** m/2n*l/21 

Si “ ' nv 17n T(l/2a) 


(A. 11) 


The problem that remains in determining the intensity is to express 
S(P) in terms of x. The rays from a source at the surface propagate, 
initially, vertically; therefore, for a vertical -force source, most of 
the energy will be directed along the z-axis with little or no energy 
propagating along the surface ■. The rays received at the surface near the 
source correspond to large values of the ray parameter p because 

p = 1/V(z t ) (A* 12) 

Where Zj is the turning depth of the ray; i.e., its maximum depth of 
penetration. Therefore, if We assume an asymptotic expansion for S(p) 
of the form (for p»i) 

S(P) ~ p" m (l + i/p + 1/p 2 + ...); m>l (A. 13) 

and use the fact that (from equation A. id) 
p *x" n 


we have 

S(p) - x 1 ™ 1 . (A. 14) 

This indicates that the source radiation pattern contributes 
an increase to the intensity as the separation increases. This is due 
to the fact that the rays detected at the larger separations crane from 


30 



that. part .of the radiation pattern vihere the intensity is higher; namely,, 
closer, to the t-axis*. In Order to insure integrability of S(p) for the_ 
large values of p (l<<p<«) , the exponent m must be greater than 1- 
Using the fact that 

1 d ! t - (n+2) (A*15) 

* 

the intensity variation with separation x becomes 


I(x) - x ‘t2*n(i-»)] 
and the amplitude variation is 

A(x) - I 1/2 - 

For the self-cdmpacting-powder model, n=l/6 aiid we have 
A(x) = A(x q ) (x/x 0 )’^ 3 " m ^ 12 ; m>l. 


(A. 16) 
(Ai 17) 
(A<18) 



Appendix B 

the relative geophdne sensitivities * shot strengths end the amplitude 
variation; least-squares analysis. 

The optimum values, in a least-squares sense, of 1) the relative geophone 

sensitivities, 2) the relative shot strengths and 3) the exponent of the 

amplitude variation with distance can be determined if there is sufficient 

redundancy in the amplitude data and the functional form of the amplitude 

function is known. The direct -wave-amplitude variation has the functional 

form given by equation (6) both for a half space (i.e. , constant velocity) 

and for a vertically inhomogeneous medium (i.e. , v(x,y,z) = v(z) only) . 

2 

The summed, weighted and squared error, E , between the log of the 
measured amplitude values and the values predicted by the functional form (as 
expressed in equation (8)) is (see equation (9)) : 

H 2 (a 0 ,m,g,s) - l l w i3 (a d + mX i;j + gi + s.. -a.,.) 2 (B.l) 

i*l j^l 

where a Q =log(GjSj) , Gj is the reference-geophone sensitivity, Sj is the 
reference-shot strength, g is the (vector of) relative geophone sensitivities, 
gi is the relative sensitivity of the i-th geophone, s is the (vector of) 
relative shot strengths, Sj is the relative strength of the j-th shot, is 
the log of the separation between the i-th geophone and the j-th shot and 
a _ is the log Of the amplitude measured at the i-th geophone for the j-th shot. 

We define a parameter vector, p, in the parameter Space made up of 
the a 0 , g^ s. and m: 
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I+J+2 


P-I g A ♦- 1 SjSj + m ♦ Vo • I P n ?n- St5 * S * 3 o (B,2) 

i -1 j*l rt-1 

where the symbols with Carets above them-Ce.g: . ig) «• “ lt vectofs ln the 
parameter space. They are jssumed to be orthogonal Cor independent), that Is: 

k-L ■ <W * ii-*u-%> V»‘V 'fc'V 


m-a 0 " V®1 - V*j * g i -S j ' °' 


(Bo 3) 


where , for example , 6^ is 


the Kronecker delta which is defined as <$^=0 


A A 


(B.4) 


if ij*k, 6 ik «l for i=k. Therefore, 

a 0 ♦ ♦ «i ♦ s P-(* 0 + "*ij + ®i + S 3 )4 

The summed, weighted and squared error, E 2 (p), will have a minimum in the 
parameter space where its gradient is zero; that is: 

r*2 


gi ♦ • V C *° * ♦ ii * *l> ■ 0 

(B.5) 


f l w i j ( a o + + 


or 


A A 


x x Vij(^*«*i3 + 8 i + V 

i 3 


(B.6) 


- ■$> 


i J 


Where the term Inside the brockets Is a square, symmetric maxtrlx obtained 
by the dyadic product of the two vectors (See equation (B. 7), next page). The 

resulting set of equations can be written in matrix form as 

_ “ (B. 8) 

a * A*p 




where a is the vector 


5 ‘ n * nV 

■ I | 4 * "*ij * *i + s j' CB ' 9) 

the square, symmetric matrix A is defined in equation (B.7) and the vector 
p is defined in equation (B.2). 

the parameter values (a 6i m,g i and s^) are just the components of the 
vector p and they are determined by inverting equation (B.8) to give 

pa A' 1 *a. CB.10) 

This is the solution to the problem of determining the optimum (in a 
least-squares sense) parameters and it can be seen from equation (B.iO) that 
the accuracy of the solution depends upon the stability of the inverse of 
the matrix A and the errors in the vector a. these, in turn, depend upon the 
accuracy of the measurements of the separations, X^, and the amplitudes , a^, 
and the values of the weights, w^. The weights themselves are established 
by the accuracies of a^ and X^. In the analysis of the amplitude data from 
the Apollo- 16 ASE, it was assumed that the separations were measured with 
high accuracy; consequently, the weights were established only on the basis of 
the accuracy of the amplitude measurements. 

From the foim of matrix A (equation (B.7)), it can be seen that it will 

have a stable inverse if there is high redundancy in the data; that is, the 

suns of the weights over the gsophones (subscript i) and over the shot strengths 

(subscript j) as well as the double sums (over i and j) have large values. This 

will occur when the amplitude from each shot is measured accurately at each 

geophone i that is, all the w^-1. Unfortunately, this is not the case in the 
Ap011o-14 or ApollO- 16 ASE’s. 
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TABLE I 

Correspondence of binary data values (B.D.) 
with the geophone voltage (V) 
(Geophone 1, Apollo 16) 


V 

B.D. 

V 

B.D. 

V 

299 

11 

-.00363 

22 

+.02101 

279 

12 

-.00202 

23 

+.03783 

7115 

13 

-.00112 

24 

+.06813 

3958 

14 

-.00047 

25 

+.1227 

2202 

15 

-.00000 

26 

+.2209 

,1225 

16 

+.00048 

27 

+.3978 

,06817 

17 

+.00111 

28 

+.7164 

,03793 

18 

+.O02OO 

29 

1.290 

,02110 

19 

+.00360 

30 

2.323 

.001174 

20 

+.00648 

31 

4.183 


.000653 


21 


+.01167 



table LL 

Shot-to-Geophone Separations 


Separation 

No. Of 

Apollo- 14 Shot Nos.* 

Apollo- 16 Shot Nos.* 

(ft) 

(m) 

traces 

GP-1 

GP-2 

GP-3 

GP-1 

GP-2 

GP-3 

0 

0.00 

6 

21 

11 

1 

19 

11 

1 

15 

4.57 

5 

20 

12 

2 

- 

10 

2 

30 

9il4 

7 

19 

i3 

3 

18 

9,12 

3 

45 

13.71 

6 

18 

- 

4 

17 

8,13 

4 

60 

18.29 

6 

17 

7 

- 

16 

7,14 

5 

75 

22.86 

4 

- 

- 

- 

15 

6*15 

6 

90 

27.43 

6 

- 

17 

7 

14 

5,16 

7 

105 

32.00 

6 

- 

4,18 

- 

13 

4,17 

8 

120 

36.58 

7 

13 

3*19 

- 

12 

3,18 

9 


*Thumper-shot numbers which had the proper separation from the 


three geophones. 
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Table Hi 

Traveltimes (ndlliseconds) 


oration Measured Traveitimes^ Calculated Traveltimes 

x(ift) 12 3 4 5 A B C 


4.57 

55 

53, 

56 

52 

- 

51.7 

44*0 

49.1 

9.14 

91 

91 

- 

87 

99 

92.1 

8ZJ 

80.2 

13.71 

123 

123 

124 

- 

128 

129.1 

131.9 

120.3- 

18.29 

151 

149 

152 

- 

155 

164.2 

175.8 

160.4 

22.86 

- 

- 

- 

- 

177? 

197.7 

219.8 

200.5 

27.43 

206? 

230 

i96? 

229 

199? 

230.1 

*245.0 

240.6 

32.06 

255? 

274? 

264? 

274? 

- 

261.7 

*260.3 

280.7 

1-n 

.76 - 

-86 

.74 

-84 

.63 

5/6_ 

— 

1 

v„(fti/s) 

590 

430 

630 

340 

1200 

350 

104 

114 


§ Times with question marks (?) indicate difficult time determinations 
* Traveltimes of the first refracted wave (earliest arrival) 

1 j Measured from Apollo- 14 and -16 stacked data (3-66 Hz) 

2. Measured from Apollo- 14 (oniy) stacked data (3-66 Hz) 

3. Measured froirt Apollo-16 (only) stacked data (3-66 Hz) 

4. Measured from Apollo- 14, geophone-2 profile (3-66 Hz) 

5. M eas ured from Apollo-16, geophone-1 profile (3-66 Hz) 

A. Self-cOnpdctihg-powder model; t«14.57 (msec.) 

B. Apollo- 14 layered model (Kovach and Watkins, 1973) 

C. Apolio-16 layered model (Kovach and Watkins, 1973) 
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Figure 1. 
Figure 2. 
Figure 3. 

Figure 4. 

Figure 5. 
Figure 6. 

Figure 7. 

Figure 8. 

Figure 9i 
Figure 10. 
Figure 11. 


List of Figure Captions 


Long-period, vertical component (LPZ) lunar seismograms, 
Apollo- U Passive Seismic Experiment, 1969. (From Latham, 
et.al-*. 1970a). 

Seismic ray paths for a linearly increasing velocity variation 
with depth and topographic irregularities. (From Cold and 
SOter, 1970) » 

Plan, view of the geophone siting and thumper-shot locations 
for the Apollo- 14 and Apollo- 16 Active Seismic Experiments* 
(Geophones : o; shots: x; misfired shots: »). 

Representative raw-data traces from the Apollo- 16 ASE. 

(Thumper shot 10). 

' 'Degl itched*vers ions of the traces in Figure 4* 

Bandpass - filtered versions Of the traces in Figure 4. (3db 
frequencies! 10 .£ § 66.25Hz). 

Amplitude spectrum of the first two seconds Of the signal 
from geophone 2, thumper shot Id, Apollo- 16 ASE. (Separation: 
4.57m). 

Stacked, filtered and amplified traces for shot-to geophone 
separations of 0, 4.57 and 9. 14nt. (Apollo-14 and Apollo- 16 
signals combined; Bandpass: 10.5 to 66Hz). 

Stacked, filtered and amplified ASE profile (Apollo- 14 and 
Apolid*16 signals combined; bandpass : 20 to 50Hz). 

Single- geOphone profile, filtered and amplified. (Geophone 1, 
Apollo- 16 ASE; Bandpass: 20 to 50Hz) . 

Log- log plot Of the traveltimes versus separations for the 
stacked and filtered traces. (Apollo- 14 and Apollo- 16 ASE 
Signals combined; Bandpass: 3 to 66Hz; measured slope: m*0.76 
and reference velocity: v Q 6 590m/sec) . 
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lUgure 12. 


Figure A-l. 


List- of Captions (Can't.} 

Log- log plot of the amplitudes versus separations (Single 
geophonc amplitudes; Geophone. 3, Apollo- 16 ASE; figndp ass: _ 

3 to 66Hz l measured slope : -2. 01} . 

a) Ray paths for a power- law velocity, variation: v(z)«v 6 (z/z 0 ) . 

b) Detail of the wavefront, dw, iri a ray bundle. 


LM impact 


DeWmber 10 


December 10 


Time, min 


Figure l. Long-period, vertical component (LPZ) lunar 
seismograms Apollo* 11 Passive Seismic Experiment ^1969. 
(From Latham, et.al. * 1970a). 




TO SEISMOMETER 




guTe 2. Seismic ray paths for a linearly increasing 
‘loeity variation With depth and topographic irregularities 
•'rom Gold and Soter, 1970). 




SHOT NO. 

2.1 20 Id 18 17 16 15 14 13 12 II 10 9 8 7 6 8 4 3 2 I 

X X X X X (g) ® ® X X X ® ® ® X ® ® X X X X 

O H-H O o 

eop I 4.57m GeOp 2 Geop 

H 45.72 m 45.72 m H 

a) APOLLO -14 ACTIVE SEISMIC EXPERIMENT. (21 shots) 
Shots 5, 6, 8, 9, 10,14,15 and IS misfired. 


SHOT NO. 


19 

18 17 16 15 14 13 12 

11 10 9 

8 7 6 5 4 3 

2 1 

X 

X X X X XXX 

X X X 

X X X X XX 

X X 

0 

H-H 

o 


o 

Geop 1 

4 . 57m 

Geop 2 


Geop 

L^ . 

45.72 m 



45.72 m 

— M 


b) APOLLO -16 ACTIVE SEISMIC EXPERIMENT. (19 shots) 


Figure 3. Plan view of the geophone siting and thumper- shot locations for 
the Apollo- 14 and ApollO s 16 Active Seismic Experiments. (Geophones: o; 
shots: x; misfires: A. 







Figure 6. Bandpass- filtered versions of the truces in Figure 4. 
(3 db frequencies: 10.5 $ 66.25 Hz). 



AMPLITUDE SPECTRUM: 



figure 1 . Amplitude spectrum of the first two seconds of the 
signal from geophone 2, thumper shot 10, Apollo-16 ASE. 
(Separation i 4.57 m). 
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STROKED, FILTERED AND AMPLIFIED 




Figure 8. Stacked, filtered and amplified traces for Shot-to 
geophone separations of 0, 4.57 and 9.14 m. (Apollo- 14 and 
Apollo- 16 signals combined; Bandpass: 10.5 to 66 Hz). 




STACKED, f IuTEAEO AND RMPLIfltfl <20H|-S0Hi) 




Figure 9. Stacked, filtered and amplified ASH profile (Apollo- 
and Apollo- 16 signals combined i Bandpass: 20 to 50 Hz). 
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I i«ure 12. Log-log plot of the amplitudes versus separations (Single 

geophone amplitudes; (teophone 5, Apollo- 16 Atfli; Bandpass: .1 to 06 II- 
measured slope: - 2 . 01 ). 






